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Abstract — This paper presents hybrid integration between 
aerodynamic, airline procedures and other BADA-based (Base of 
Aircraft Data) coefficients with a classical aircraft dynamic 
model. This paper also describes a three-dimensional guidance 
algorithm implemented in order to produce commands for the 
aircraft to follow a flight plan. The software chosen for this work 
is MATLAB®. 
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I. Introduction 

Unmanned Aircraft Systems (UAS) have played an 
important role in the global Aeronautics industry for the past 
few decades. Even so, most UAS operations within the 
National Airspace System (NAS) are restricted due to the 
possible risk that grand scale civil UAS operations may pose 
to the safety of aerial operations. There are currently few 
military UAS operating in the NAS. Without an onboard pilot, 
UAS lose the ability to “see and avoid” other aircraft, it is for 
this reason that having accurate models is necessary. 

There are currently several trajectory generators to support 
Air Traffic Management (ATM) applications and research. 
There are also a couple similar to what is here presented, such 
model as KTG [6] and HYBRIDGE [2]. However, this model 
has as objective developing a stand-alone MATLAB tool in 
which the performance and behavior of several guidance, 
collision avoidance and sense and avoid algorithms can be 
observed. This will aid research in Next Generation ATM. 

II. BADA Models 

BADA provides a set of ASCII (American Standard Code 
for Information Interchange) files that contain performance 
parameters and coefficients for 399 types of aircraft. This data 
is used by EUROCONTROL to simulate standard operations 
on European airspace. There are present efforts in creating 
BADA compatible files for several unmanned systems. A 
brief explanation of the most important BADA models and 
parameters is now given. Refer to [1] for a broader 
explanation on each model. 

A. Atmosphere 


This model provides expressions for temperature, pressure, 
density and speed of sound as a function of altitude. These 
expressions are based in the International Standard 
Atmosphere. It is dived by the tropopause, defined at 11000 m 
in terms of geopotential pressure altitude. This model also has 
temperature and pressure corrections for non-standard 
atmospheres. 


B. Total Energy Method 


Total Energy Method relates de rate of work done by 
forces acting on the aircraft to the rate of change in potential 
and kinetic energy. It is given by the following expression: 
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Where Thr means thrust, D is the aerodynamic drag force, 
V TA s is True air speed, m means the mass of the aircraft, g 0 is 
gravitational acceleration and finally h is geodetic altitude. 

From this equation, an expression for the Rate of Climb or 
Descent (ROCD) can be obtained. 


C. Aerodynamics 

BADA contains a simple Aerodynamic model in which the 
aircraft’s angle of attack is always considered to be zero; there 
is a correction for bank angle nevertheless. Drag coefficient 
depends on flight configuration. Cruise drag coefficient is 
valid for all flight configurations except approach and landing, 
where additional coefficients must be used to simulate the use 
of flaps and extension of landing gear. 


D. Thrust Model 

Thrust produced by the aircraft depends directly on the 
engine type, which can be jet, turboprop or piston. 
Additionally, flight configuration such as Climb (CL), Cruise 
(CR) or Descent (DES) has a direct impact on the amount of 
thrust available. Several other parameters like altitude and 
speed influence this thrust model. 

E. Fuel Consumption 

Fuel consumption is a direct result of the aircraft’s thrust, 
which implies a direct relationship with engine type and flight 



configuration. Thrust Specific Fuel Consumption (TSFC) is 
determined via engine coefficients, speed and altitude. 

F. Airline Procedure Models 

Airline procedure models contain the speed schedules for 
CL, CR, and DES configurations and dictates Calibrated Air 
Speed at different altitudes. For each of these phases and each 
aircraft model, Airline procedure model requires default 
speeds, stall speed for take off and landing configurations and 
several other coefficients such as maximum operating Mach 
number. Figure 1 shows the speed schedule for a CR 
configuration of an Airbus A3 00 aircraft. 


Speed Schedule 



Fig.l Speed Schedule on a cruise configuration for an Airbus A300 aircraft 


G. Global Aircraft Parameters 

BADA incorporates many parameters that are independent 
of aircraft type known as Global Aircraft Parameters. For 
instance, standard operational values such as maximum 
normal and longitudinal acceleration, nominal and maximum 
bank angle for several flight configuration, minimum speed 
coefficients and thrust coefficients are established in Global 
Aircraft Parameters. 

III. Aircraft’s Equations of Motion 

The aircraft dynamic model is given by three control inputs 
and seven differential equations, which represent the aircraft’s 
states. The state vector is given by: 

X (?) = [x(?), y(t\ h(t),V TAS (?), i//(/), m{t), < o{t )] 

Table 1 presents the system’s control inputs as well as the 
equations of motion. 


TABLE I 

Control Inputs and Equations of Motion 
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T and (j) command are the guidance algorithm outputs, ys is the 
aircraft’s heading with respect to the north, the temperature 
factor ( T-AT/T) is a correction for non standard atmospheres, 
f{M} is known as the energy share factor and it specifies how 
much of the available power is allocated to climb as opposed 
to acceleration while following a selected speed profile during 
climb, it has been shown that it is a function of Mach number, 
(refer to [6] for a broader explanation), (j) is the aircraft’s bank 
angle, rj is the fuel consumption and k is the bank gain. 

Thrust in this case is more of a system constraint than a 
control input; there is no direct throttle input. The equation for 
ROCD is obtained from the TEM presented in previous 
sections. 

Aircraft’s bank response is approximated as a first order 
response, while flight path angle is considered immediate. 
Future work on this model may improve these responses. 

Change in heading is assumed as level-turn flight. Change 
in mass in the aircraft is only given by the fuel consumption. 

Additional to these aircraft states, the variables 
“climbmode” and “accelmode” were created; these 
variables take the values of CL, CR and DES, and accelerate 
(A), cero acceleration (C) and deceleration (D) respectively. 
These values indicate the flight configuration and acceleration 
mode required for the aircraft to follow the predetermined 
speed schedule established in BADA’s Airline Procedure 
Model. These variables have a direct impact on thrust, fuel 
consumption, aircraft speed and ROCD. 

IV. 3D Guidance Algorithm 

This algorithm outputs the necessary commands for the 
aircraft so that it follows a predetermined flight plan given by 
waypoints. The algorithm uses an X-Y-Z frame of reference in 
which positive X direction represents east, positive Y 
represents north and positive Z represents altitude. Guidance 
logic is non linear and creates lateral acceleration commands 
that can be easily translated into bank angle commands. This 
logic is based on the work presented in [3], the waypoint 
tracking and waypoint target switching presented in [4] is also 
used. 

In order to have a complete 3D guidance algorithm, a flight 
path angle command was added in the vertical plane. Its 
derivation is purely geometrical and is simply expressed as: 


y = tan 

Where zl/z is the altitude difference between the aircraft’s 
current position and the altitude of the target waypoint and 
Aground represents the distance on the horizontal plane between 
the aircraft and the target waypoint. 

Maximum normal acceleration values established in 
B ADA’s Aircraft Global Parameters must be taken into 
account. This value will limit the changes in flight path angle 
command and thus rate of climb or descent. This serves the 
purpose of simulating realistic aircraft behavior according to 
standard operations. 

V. B ADA, Equations of Motion and Guidance 
Algorithm Integration 


Ah 
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was compared to 10 flights in ACES background traffic to 
verify it performs adequately. Figures 3 and 4 show the RMS 
error of the displacement in both the vertical and horizontal 
planes. 



This section contains a brief explanation on how the 
different models are integrated forming the hybrid model. 

Figure 2 shows a block diagram indicating the flow of 
information within the models. The process begins with a pre- 
established flight plan that consists of a series of waypoints, 
the guidance algorithm takes this information as well as the 
aircraft’s initial position to determine a target waypoint and 
compute the necessary control commands. 



Fig. 2 Block Diagram of the Hybrid Model indicating the flow of information 

The block labeled as “BADA Models” takes information 
from the initial condition and the control commands computed 
by the guidance algorithm to determine values of speed 
schedule, thrust, drag, fuel consumption, climb mode, 
acceleration mode and others. All this information is fed into 
the block labeled as “Equations of motion” which, as the 
name implies, contains the equations of motion that are solved 
using the Euler method with a time step of 1 second. The 
updated aircraft states are then fed to both the guidance 
algorithm and BADA models to determine the next set of 
commands. The hybrid trajectory generator here developed 



It is worth mentioning that the mean distance between 
waypoints in these comparisons was approximately 13.23 
nautical miles. It can be seen that the RMS error in the vertical 
plane is much greater than the error in the horizontal plane. 
This can be explained by the difference in the way ACES 
solves the climb and descent phases compared to the model 
here presented. While ACES takes in an altitude clearance and 
solves a boundary value problem, the hybrid trajectory 
generator takes in specific altitudes for each waypoint. 
Nevertheless, the mean RMS error in the horizontal plane for 
the 10 flights is approximately 17.6 ft. This is an acceptable 
result. 

VI. Conclusions 

This document showed a hybrid trajectory generator, which 
combines models, parameters and coefficients from the Base 
of Aircraft Data and a classic aircraft dynamics. Also, the 
implementation a 3D guidance algorithm was presented. 

This model will be a useful tool to simulate the interaction 
between conflicting aircraft according to aircraft’s 

performance as well as standard operations. Present and future 
UAS research may benefit from this model. It has applications 
in Sense and Avoid research as well as testing the 
performance of new control and guidance algorithms. 
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